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The large-scale circulation of planetary atmospheres like that of the Earth is traditionally thought 
of in a dynamical framework. Here, we apply the statistical mechanics theory of turbulent flows 
to a simplified model of the global atmosphere, the quasi- geostrophic model, leading to non-trivial 
equilibria. Depending on a few global parameters, the structure of the flow may be either a solid- 
body rotation (zonal flow) or a dipole. A second order phase transition occurs between these two 
phases, with associated spontaneous symmetry-breaking in the dipole phase. This model allows us 
to go beyond the general theory of marginal ensemble equivalence through the notion of Goldstone 
modes. 
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I. INTRODUCTION 

It is well known that two-dimensional (2D) turbu- 
lent flows develop large-scale coherent structures [T] , and 
widely believed that the existence of these structures is 
related to the existence of additional conserved quanti- 
ties, as compared to the 3D case. In particular, Kraich- 
nan [5] showed that in 2D flows, due to the conservation 
of enstrophy and energy, two inertial ranges coexist: a 
direct (downscale) enstrophy cascade and an inverse (up- 
scale) energy cascade. As a result, the energy accumu- 
lates at the largest scales (lowest wavenumbers) leading 
to coherent structures. This phenomenon is sometimes 
related to a form of Bose-Einstein condensation [21 [3]. 
The formation of these coherent structures can be ac- 
counted for by statistical mechanics approaches. Given 
some conserved quantities, one is able to predict the most 
probable flow, which, assuming ergodicity, coincides with 
the final state of the flow. One particular area where 
large-scale coherent structures are ubiquitous, and dy- 
namical approaches are both analytically intractable and 
numerically not affordable, is geophysical fluid dynamics. 
The early attempts to gain qualitative understanding of 
geophysical flows from statistical mechanics focused on 
the spectral space jH [5] , but a theory dealing with con- 
tinuous fields in the physical space is now available [SHE] . 
The Millcr-Robcrt-Sommeria (MRS) theory was used to 
develop small-scale parameterizations [51TH] and study 
the formation of localized vortices like Jupiter's great red 
spot [151417] or meso-scale ocean structures like vortex 
rings and zonal jets [T51 [TS]. It has also been discussed 
in relation to Fofonoff flows [2"DH2"2] . 

In this paper, we construct the phase diagram of a 
simple model of atmospheric flow on a rotating sphere, 
the barotropic vorticity equation. Non-trivial equilibria 
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are obtained and a second-order phase transition with 
spontaneous symmetry breaking is observed between a 
solid-body rotation phase and a dipole phase. In ad- 
dition, the system displays surprising thermodynamic 
properties. Standard thermodynamics textbooks usually 
explain that for large enough systems (thermodynamic 
limit), mean values computed with a probability distri- 
bution that assigns equal weights to all the states hav- 
ing a given energy (microcanonical distribution) coincide 
with mean values computed with a Boltzmann distribu- 
tion (canonical distribution). When the energy is not 
additive, which occurs in the presence of long-range in- 
teractions [23l [24] , this is generally not the case. Con- 
sequently, peculiar thermodynamic features may appear: 
for instance negative microcanonical specific heat [551 I2S] 
(due to the fluctuation-dissipation theorem, the canonical 
specific heat is always positive), convex dips in the mi- 
crocanonical entropy [27] and unusual phase transitions 
[2"51 12"5] . The statistical equilibria of the barotropic vor- 
ticity equations on the rotating sphere have been the ob- 
ject of a number of studies in the past [30U321 . However, 
no complete theory in the general framework of modern 
statistical mechanics of turbulent flows is available. The 
purpose of this paper is to fill this gap by providing the 
phase diagrams in the various statistical ensembles and 
discussing in detail a notion of marginal ensemble equiv- 
alence. 



II. STATISTICAL MECHANICS OF IDEAL 
FLOWS ON A ROTATING SPHERE 

We base our study on the quasi-geostrophic (QG) 
equations which provide a simple model of the atmo- 
spheric circulation. For simplicity, we shall treat here 
only the barotropic case with no bottom topography so 
that the potential vorticity is just the sum of the relative 
vorticity and the planetary vorticity: q = lu + /, where 
/ = 2fi cos and is the angular velocity of the Earth. 
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The general case will be described in detail elsewhere [53] . 
In the absence of forcing and dissipation, the potential 
vorticity is advected by the flow: 

d t q + u ■ Vg = 0. (1) 

Here u = — f x Vtp is the velocity field and ip the stream 
function, related to the vorticity field by oj — —At/;. For 
convenience, we adopt the Gauge condition J^ip dr = 0. 
Similarly to the 2D Euler equations, QG flows admit dy- 
namical invariants: the energy E[q] — (1/2) f v {q— f)ip dr 
and the Casimirs T n [q] = J v q n dr. For specific flow 
domains V, some additional invariants must be con- 
sidered. Here, we focus on the case of a sphere with 
unit radius T> = S 2 , for which the angular momentum 
L[<l] = Jr>(<7 — /) cosOdr is conserved. 

Steady-states of the QG equations which are stationary 
in a frame rotating with angular velocity Ql with respect 
to the terrestrial frame satisfy q(9, <j>, t) = q(9, <f> — fi.Lt). 
Substituting this relation into the evolution equation (HI , 
we obtain {q,ip + QlcosB} = where {,} denote the 
Poisson brackets. Hence the general form of steady states 
of the QG equations on the sphere is q — F{ip + £l L cos 9) 
where F is an arbitrary function. Statistical mechanics 
allows us to select the most probable function F given the 
value of the dynamical invariants E, L and T n as follows. 

A priori, the evolution of the flow is purely deter- 
mined by the initial potential vorticity field q$. It is well 
known that throughout the dynamical evolution, the po- 
tential vorticity field develops finer and finer filaments. 
This small-scale behaviour of the potential vorticity field 
makes it difficult to treat the flow in a deterministic way, 
as the typical scale of potential vorticity fluctuations will 
eventually become smaller than any simulation or ob- 
servation resolution. On the other hand, the large-scale 
structure can be readily predicted by statistical mechan- 
ics. Introducing a coarse-grained potential vorticity field, 
one maximizes the statistical entropy subject to the con- 
servation constraints to obtain the most probable poten- 
tial vorticity field [SHE] . It is characterized by a functional 
relationship q = F(ijj + Hlcos9) where the function F 
depends on the invariants. In the following, we shall as- 
sume that F is linear. Several justifications can be given 
to motivate this choice: (i) it corresponds to a MRS equi- 
librium with a Gaussian potential vorticity distribution 
reached for particular initial conditions 6J. Furthermore, 
a state that maximizes S[q] = —(1/2) J v q 2 dr at fixed 
energy, circulation and angular momentum is granted to 
be thermodynamically stable [20l|34] (see also [H])i (ii) a 
linear q — ip relationship is obtained in the limit of strong 
mixing and in the low energy limit |35| : (iii) maximizing 
the statistical entropy with conservation of energy, cir- 
culation, angular momentum and fine-grained enstrophy 
is equivalent to minimizing the coarse-grained enstrophy 
(1/2) J v q 2 dr defined in terms of the coarse-grained po- 
tential vorticity at fixed energy, circulation and angular 
momentum [36] ■ The ensemble-mean potential vorticity 
is characterized by a linear q— ip relationship and the fluc- 
tuations around it are Gaussian; (iv) some authors have 



proposed to treat the Casimirs in a canonical way [37], 
which is equivalent to choosing a prior small-scale poten- 
tial vorticity distribution [57H35] . The ensemble-mean 
potential vorticity then maximizes a generalized entropy 
S[q] = — J v C(q) dr where C is a convex function related 
to the prior. A Gaussian prior leads to a linear q — tp 
relationship; it is associated with a quadratic general- 
ized entropy S[q] = —(1/2) J v q 2 dr proportional to the 
coarse-grained potential enstrophy. All these justifica- 
tions lead to a linear q — ip relationship. However, they 
are not equivalent. They essentially differ in the manner 
that the fragile constraints (high-order Casimir invari- 
ants) are treated, while the robust constraints (energy, 
circulation, angular momentum) are strictly taken into 
account in all these approaches. In approaches (i) and 
(ii), although the structure of the equilibria in the limit 
considered depends only on the robust invariants, con- 
servation of the fragile invariants is taken into account 
implicitly. On the contrary, approaches (iii) and (iv) are 
based on the idea that, in real flows, forcing and dissipa- 
tion acting at small scales break down the conservation of 
high-order Casimirs. In (iii) these constraints are purely 
discarded while in (iv) they are replaced by a prior. Note 
that strictly speaking, either we consider the unforced- 
undamped case and we should consider all the invariants, 
or we consider the forced-damped case in which case the 
value of the robust invariants should also be selected by 
the forcing and the dissipation. Finally, the approach (ii) 
does not make any assumption on the fluctuations of the 
potential vorticity, contrary to (i), (iii) and (iv) in which 
they are Gaussian. 

The microcanonical variational problem reads 

S(E, L) = m^{S{q}\E[q] = E, L[q] = L}. (2) 

9 

Note that it is not necessary to include the circulation 
in the constraints as it always vanishes due to the ge- 
ometry. The solutions of this variational problem are a 
subset of the full class of MRS equilibria. In practice, 
constrained variational problems like ^ are difficult to 
solve. It is much more convenient to work with the dual 
variational problem with relaxed constraints, the grand- 
canonical variational problem: 

J(fi, /i) = max{%] - pE[q] - M L[g]}. (3) 
i 

Clearly the two variational problems have the same criti- 
cal points, but the nature of the critical points (maxima, 
minima, saddle points) may differ. When this happens, 
we speak of ensemble inequivalence |27l 1371 140] . For iso- 
lated systems, like the large systems encountered in as- 
trophysics or geophysical fluid dynamics, the natural en- 
semble is the microcanonical ensemble. The classical in- 
terpretation of the grand-canonical ensemble is that the 
system is in contact with a reservoir of heat and angular 
momentum. The physical interpretation of such a reser- 
voir is unclear. Nevertheless, it is always worthwhile to 
consider the grand-canonical ensemble, be it just as a 
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useful mathematical device. Indeed, a solution of the 
grand-canonical problem ^ is always a solution of the 
more constrained dual microcanonical problem ([2| . Note 
that the variational principles ^ and ^ also provide 
sufficient conditions of nonlinear dynamical stability for 
the QG equations [37] (see also [14]). 



III. SOLUTIONS OF THE MEAN-FIELD 
EQUATION AND PHASE TRANSITION 

The critical points of ^ and (|3) satisfy SS - j3SE - 
fiSL = 0, which yields a linear q — tp relationship: q = 
—/3tp — fJ- cos 9, in accordance with the general form for 
steady-states q = F(tp + Q L cos9), where L = n/P is 
the angular velocity of the frame in which the flow is 
stationary. Replacing q by its definition, we obtain the 
mean-field equation 
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FIG. 1: (color online). Amplitude of the dipole, relative to 
the solid-body rotation, as a function of the control parameter 
e. At e = e c = 1, a second-order phase transition occurs, ac- 
companied by spontaneous symmetry breaking (in the dipole 
phase). Sample stream functions are shown in the insets, with 
the arbitrary phase <po se t to 0. 



Alp - f3tp = f + n cos 9. (4) 

To solve this Helmholtz equation, we introduce the eigen- 
values and eigenvectors of the Laplacian on the sphere: 
AY nm = (3 n Y nm , where the spherical harmonics Y nm 
form an orthonormal basis of the Hilbert space L 2 (T>), 
and f3 n = —n(n + 1) (n G N, — n < m < n). 

When f3 is not an eigenvalue of the Laplacian, the op- 
erator A — f3I is invertible and the equilibrium stream 
function reads 

tp = £1* cos 6, (5) 

with = (2£! + ^i)/(/3i — f3). This corresponds to a solid- 
body rotation with angular velocity O*. The energy, an- 
gular momentum and entropy depend only on the angular 
velocity: E = Q, 2 /3, L = 20*/3, S = -2(0 + 0*) 2 /3. In 
particular, for a solid-body rotation, the energy and an- 
gular momentum are not independent conserved quan- 
tities: they are related by E — E*(L) where E^(L) = 
3L 2 /4. One can easily show that the energy of any flow 
with angular momentum L is always greater than the 
energy of the solid-body rotation with the same angular 
momentum: E > E^L). 

When /3 is an eigenvalue of the Laplacian but not the 
first non-zero one - say /3 = fi n >i, the stream function be- 
longs to a 2n+l dimensional vector space: tp = fi* cos 9+ 
Z) m V'mn^nm with again = (2Q + fi)/(j3i-(3 n ). It is a 
superposition of a solid-body rotation with a multipole. 
The energy, angular momentum and entropy are found 
to be E = - (3 n J2 m Hnra\ 2 /2, L = 2Q*/3, S = 

-2(0 + O073 " Pn Em IVw| 2 /2- 

Finally, the mean-field equation admits solutions for 
P = fix only if fi = fi c = —20. In that case, the stream 
function reads tp = 0„ cos 9 + 7 C sin 9 cos <p + 7s sin 9 sin (p. 
The energy, angular momentum and entropy depend on 
the coefficients 0*,7 C ,7 S : E — (O* 2 + 7 2 + 7 2 )/3,L = 
20*/3, S = -2((0 + 0,) 2 +7 2 + 7 2 )/3. Introducing the 
angle (po such that 7 C = y/3(E — E*(L)) cos (po and 7 S = 



\/3(E — E*(L)) sin 0o, the stream function becomes 
tp = Q„ cos 6 + v/3(£ - sin 6 cos(0 - <p ). (6) 

This is a dipole flow with an arbitrary phase (pa. Setting 
e = E/E*(L), the amplitude of the dipole relative to the 
background solid-body rotation can be recast as a(e) = 
\Jt — 1; then tp = (3L/2)[cos6> + a(e) sin 9 cos((p — (po)]- 
This reveals a second-order phase transition between a 
solid-body rotation phase at e = e c = 1 and a dipole 
phase for e > e c (Fig. [T]). The angular momentum con- 
trols the overall amplitude. The phase cpy is not deter- 
mined by the control parameters: this is a case of spon- 
taneous symmetry breaking. Since the phase transition 
line coincides with the line of minimum energy E„(L), 
the phase transition would be difficult to see in practice 
(in numerical simulations for instance). This is reminis- 
cent of some systems encountered in condensed matter 
physics, where the critical point is sometimes reached at 
vanishing temperature (T = 0), like for instance in the 
ID Ising chain. However, it is possible that, in more 
realistic models, the two lines separate from each other. 

We stress that the spontaneous symmetry-breaking 
property is also related to an apparent paradox of the 
Kraichnan statistical theory in spectrally truncated space 
|41) . In the spectrally truncated statistical mechanics of 
quasi-geostrophic flows over topography [?J [SJ S2] , there 
is usually a distinction between the mean flow and the 
eddies, or transients. Statistical equilibria are said to be 
sharp when the eddy component vanishes in the limit of 
infinite cut-off wavenumber. An often-read claim is that 
in the absence of topography, all the energy is contained 
in the transients and there is no mean flow in the sta- 
tistical equilibria, breaking the sharpness property. This 
would be in contradiction with the MRS theory which 
directly predicts sharp equilibria. This apparent para- 
dox is waived by a symmetry breaking argument: for a 
given value of the canonical parameters, opposite stream 
functions are equally admissible. When summing over 
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all the states, they cancel out, leaving a null mean flow. 
However in practice, the system selects spontaneously a 
particular symmetry-breaking state within the symme- 
try restoring set of equilibrium states. Note that this 
apparent paradox is more easily understood in the con- 
text of inequivalence of statistical ensembles: the proper 
ensemble to work with is the microcanonical ensemble, 
in which the paradox does not appear. But this was not 
noticed at that time because statistical mechanics in the 
microcanonical ensemble was still poorly understood (the 
first observation of ensemble inequivalence was made in 
[251 I2D] , but it took some time before a thorough under- 
standing was reached [27]), although Kraichnan himself 
expressed concerns about the relevance of canonical mea- 
sures for fluid systems. 



IV. NATURE OF THE CRITICAL POINTS AND 
MARGINAL ENSEMBLE EQUIVALENCE 

The nature of the critical points of entropy is given 
by the second variations of the grand-potential func- 
tional J[q] = S[q] — 0E[q] — ftL[q]: a critical point 
of entropy at fixed energy and angular momentum is 
a local maximum if and only if S 2 J < for all per- 
turbations 5q that conserve E and L at first order 
[14j . This is the microcanonical stability condition. In 
the grand-canonical ensemble, the stability condition is 
S 2 J < for all perturbations. The quadratic form 
8 2 J = -(l/2)J v (5q) 2 dr- (1/2)0 J v {\75^) 2 dr is neg- 
ative definite for > 0%. Hence, in this case, the flow is 
grand-canonically stable and therefore also microcanoni- 
cally stable. When < 0\, it is easy to build perturba- 
tions that destabilize the flow. Consider a perturbation 
8q proportional to a Laplacian eigenvector Y nm . The ba- 
sic flow is either a solid-body rotation, or a degenerate 
flow 03 = P>1 ). For (n,m) ^ (0,0), (1,0) (and n ^ p 
in the degenerate case), the perturbation conserves the 
energy and the angular momentum, while at the same 
time 5 2 J > 0. Thus, all the states obtained for f3 < (3i 
are microcanonically unstable, and therefore also grand- 
canonically unstable [33] , 

When = 0\ the situation is slightly more subtle: 
the quadratic form S 2 J is degenerate: its radical R is a 
vector space of dimension 3, spanned by Yn,Yi t -i and 
Yio- All the elements of this vector space (which are 
combinations of solid-body rotations and dipoles) share 
the same value J = — 2£l 2 /3 of the grand-potential func- 
tional, which means that they are metastable in the 
grand-canonical ensemble. Perturbations belonging to 
R allow transitions between one flow of the form (|6| to 
another with different values of L,E and 4>q. In the mi- 
crocanonical ensemble, imposing conservation of angu- 
lar momentum and energy at all orders implies that the 
only perturbations in R to be considered are those which 
only change </>o- These solutions have the same entropy 
S = -2{E - £»(£)) - (3/2)(L + 2f7/3) 2 . Hence each 
state with given (E, L) in the microcanonical ensemble is 



metastable, but with fewer Goldstone bosons (only one) 
than in the grand-canonical ensemble (where there are 
three Goldstone bosons). 

To sum up, in the grand-canonical ensemble, the equi- 
librium flow is either a solid-body rotation (/3 > 0x, co- 
rotating for /J, < fi c , counter-rotating for fi > /i c ) or a 
dipole (0 = /3i, /j, — fj, c ) with arbitrary energy, angu- 
lar momentum and phase. The corresponding phase di- 
agram is represented on Fig. [2] In the particular case 
/i = (no angular momentum constraint), we recover the 
solid-body rotations obtained previously in Monte-Carlo 
simulations [3D]. In the microcanonical ensemble, when 
E = E*(L) the equilibrium flow is a solid-body rotation 
(co-rotating for L > 0, counter-rotating for L < 0), while 
for E > E*(L) it is a dipole with arbitrary phase (Fig. 
[3]). A second-order phase transition occurs between the 
solid-body-rotation phase and the dipole phase, accom- 
panied by U(l) spontaneous symmetry breaking, as the 
phase of the dipole is undetermined. The choice of a par- 
ticular phase is analogous to the choice of a fundamental 
state in field theory [43] . 

The second order phase transition is very special be- 
cause one of the two phases (solid-body rotation) exists 
only on a curve E*(L) while the other phase (dipole) 
exists on a surface E > E*(L). On the curve E*(L), 
the parameters 0(E,L) and fi(E,L) take any values in 
the range > 0i and fi < fi c (if L > 0) or fi > [i c (if 
L < 0), while on the surface E > E*(L) they take the 
unique values = 0\ and \i = /i c . Therefore, the second 
derivatives d 2 S/dE 2 = d(i/dE and d 2 S/dL 2 = d\xjdL 
pass from +00 to as we go from the solid-body rota- 
tion phase to the dipole phase. This is therefore a case 
of extreme discontinuity of the second derivatives of the 
entropy with respect to the conserved quantities. 

Note that the analytical tractability of the system 
considered here allows us to understand precisely how 
the grand-canonical equilibrium states and the micro- 
canonical equilibrium states are related: to every half- 
straight line starting from the critical point on the grand- 
canonical phase diagram (in the right half-plane, see Fig. 
[2]), with equation /i = f2*(/3i — (3) — 2il, we can associate 
a single point (E, L) = (fl 2 /3, 2f2*/3) on the E = E*(L) 
parabola in the microcanonical phase diagram (Fig. pjj). 
Similarly, the critical point (0 = pi,H = fic) is mapped 
onto the whole area over the parabola. Formally, the 
grand-canonical and microcanonical ensembles are equiv- 
alent at the macrostate level [ID]; it is in fact a case of 
marginal equivalence |27j . Note that the marginal ensem- 
ble equivalence is extreme here in that the vast majority 
of the equilibrium states in the microcanonical ensemble 
(E > E* (L)) are all included in the set of grand-canonical 
equilibria obtained at a single value of the Lagrange pa- 
rameters (0 = 0i, fj, = fj, c ). As a consequence, the con- 
struction of equilibrium states in the grand-canonical en- 
semble may be difficult to control in practice. Marginal 
ensemble equivalence is also seen at the thermodynamic 
level: geometrically, the microcanonical entropy S is a 
plane. Therefore, it is a concave function, which is char- 
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FIG. 2: (color online). Grand-canonical phase diagram 
for barotropic flow on a rotating sphere. Stable equilib- 
rium states are obtained for /3 > At the critical point 
(/3 = = /i c ), there are an infinity of metastable states. 

Vertical dashed lines (green) indicate the position of unstable 
degenerate states. At the intersection with the line n = fi c , 
(dashed red points), only the degenerate part subsists. On the 
horizontal blue solid (resp. dashed) line the equilibrium flow 
is a stable (resp. unstable) trivial motionless flow. No solu- 
tion exists if = /3i and /i 7^ /i c . The dotted half straight-line 
represents a iso-il* line. 
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FIG. 3: (color online). Microcanonical phase diagram for 
barotropic flows on a rotating sphere. The parabola E = 
E,(L) corresponds to solid-body rotations, while the shaded 
area over the parabola represents dipole flows. The area below 
the parabola is not accessible. The parabola is the location of 
a second-order phase transition with spontaneous symmetry- 
breaking. 



acteristic of ensemble equivalence at the thermodynamic 
level [23 HO], but only marginally so: it is also a convex 
function. Equivalently, the specific heats d 2 S/dE 2 and 
d 2 S/dL 2 vanish. This indicates that the system con- 
sidered is on the edge of ensemble inequivalence: the 
grand-canonical and microcanonical ensembles are for- 
mally equivalent, but very close to being inequivalent. 

Here it is possible to analyze the ensemble equivalence 
properties beyond marginal ensemble equivalence. In- 



deed, as explained above, the mctastability in the grand- 
canonical ensemble is more severe than in the micro- 
canonical ensemble since only the phase of the dipole 
is undetermined (and thus free to vary due to spon- 
taneously generated perturbations) in the microcanon- 
ical ensemble while the energy and angular momentum 
are also unconstrained in the grand-canonical ensemble. 
Hence a refinement of the notion of marginal ensemble 
equivalence can be introduced by considering the num- 
ber of Goldstone modes in each ensemble. It is equal to 
one in the microcanonical ensemble, three in the grand- 
canonical ensemble, and two in any mixed ensemble (not 
discussed). As a practical manifestation of these Gold- 
stone modes, numerical constructions of the equilibrium 
states have to be aware that spontaneous fluctuations 
may exist, the form of which depends on the statisti- 
cal ensemble considered. This property stems from the 
degeneracy of the first eigenvalue of the Laplacian on 
the sphere. Up to now, most studies have assumed 
that the first Laplacian eigenmode was not degenerate 
[20j EH HH SS], with the notable exception of pi 146]. 
Another important consequence is that in the presence 
of degeneracies, ensemble inequivalence results such as 
those observed in [2U] can collapse. 

Note also that the marginal ensemble equivalence and 
the second order phase transition observed here may not 
be robust with respect to small nonlinearities in the q — tp 
relationship. As analyzed in the case of the energy- 
circulation ensemble on a rectangular domain [45] , de- 
pending on the sign of the small nonlinearity, second 
order phase transitions can turn into first order phase 
transitions. Due to the formal similarity between the 
two systems, the situation is likely to be the same here. 



V. CONCLUSION 

In this paper, we have explicitly computed a class of 
statistical equilibria for 2D flows on a rotating sphere, 
built the corresponding phase diagrams and character- 
ized precisely the ensemble equivalence properties. Re- 
mark that all the energy condensates in the lowest (non- 
zero) order mode. It is the combination of the angu- 
lar momentum conservation and Laplacian degeneracy on 
the sphere that allows for non-trivial Bosc-Einstein con- 
densation [2J; with the angular momentum constraint, 
the energy is shared between two phases, which makes 
possible phase coexistence and phase transition. Hence 
the results presented here have a fundamentally geomet- 
ric nature. Note that in the dipole phase, the value of the 
angular momentum fixes the solid-body rotation part of 
the flow. The subtle ensemble equivalence properties of 
the system also stem from the particular geometry: the 
general inequality E > E* (L) is purely kinematic. With- 
out this inequality, ensemble equivalence would not hold. 
Nevertheless, ensemble equivalence is only marginal [27], 
and it is so in a particularly extreme way. At the ther- 
modynamic level, the entropy surface is concave, but not 
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strictly: it is locally flat at each point of the interior of the 
parabola E — E*(L). At the macrostate level, each mi- 
crocanonical equilibrium state is also a grand-canonical 
equilibrium state. More precisely, the set of macrostatcs 
Ai(E,L) obtained in the microcanonical ensemble with 
a given energy E and angular momentum L is a subset 
of the set of macrostates GC(/3, fx) obtained in the grand- 
canonical ensemble with the adequate Lagrange param- 
eters (0,n). When E = E*(L), M{E,L) = GC{fi^) 
for any (j3,fi) satisfying (20 + — 0) = 0* with 

ft, = 3L/2. For E > E*(L), M(E, L) C gC(/3 u fi c ). 
Indeed, the set £/C(/3i,Mc) contains all the dipole flows 
with any value of the energy and angular momentum in 
the interior of the parabola, while the set M. (E, L) only 
contains the dipole flows with given values of the energy 
and angular momentum. This is a loose constraint on the 
two macrostate sets. We have shown that it is possible 
to go further than the usual notion of marginal ensem- 
ble equivalence by introducing the concept of symmetry 
breaking: the symmetry breaking occurring in the grand- 
canonical ensemble is stronger than in the microcanonical 
ensemble. As usual in the theory of continuous symmetry 
breaking, this is measured by the number of Goldstone 
modes. Statistical mechanics in the microcanonical and 
grand canonical ensembles differ in their number of such 
modes (zero-energy modes in the microcanonical ensem- 
ble and zero-temperature modes in the grand canonical 
ensemble). This subtle difference can exist only when the 
Laplacian is degenerate on the flow domain, which is the 
case on the sphere and on any manifold with sufficient 
symmetries (in fact, it suffices that the eigenspace corre- 
sponding to the first non-zero eigenvalue of the Laplacian 
be degenerate). For instance, the analysis does not apply 
to the case of a rectangular domain with linear q — tp re- 
lationship and zero circulation [5TJ [35] , even though the 
entropy depends linearly on the energy in that case. Con- 
trary to our study, symmetries are not sufficient in this 
example to go further than the notion of marginal ensem- 
ble equivalence: on a rectangular domain, one only has a 
discrete, not continuous, symmetry breaking, which pre- 
vents the discussion in terms of Goldstone bosons. This 
explains why it was not discussed before, as all authors 
focused on cases when the first eigenvalue of the Lapla- 
cian is not degenerate. Therefore, there exist cases for 
which a refinement of the definition of marginal ensemble 
equivalence is needed. As Laplacian degeneracies stem 
from symmetries, it is natural that the degree of sym- 
metry breaking in macrostates becomes the refined no- 
tion of marginal ensemble equivalence at the macrostate 
level. In our system, we are able to say that there is one 
Goldstone mode in the microcanonical ensemble, three 
Goldstone modes in the grand-canonical ensemble and 
two Goldstone modes in both mixed ensembles (where 
one constraint is treated microcanonically and the other 
canonically). 

Note that in the case of a rectangular domain 
with doubly-periodic boundary condition (the two- 
dimensional torus), the first Laplacian eigenmode is de- 



generate (for aspect ratio one). The consequence is that 
the topology of the flow is not fixed by the invariants 
when considering only the maximization of a quadratic 
generalized entropy functional with energy constraint: 
both parallel flows and dipole flows belong to the fam- 
ily of extrema [46] . The degeneracy can be lifted either 
by stretching the domain or by considering small non- 
linearities in the q — ip relationship. The coexistence of 
two phases is very similar to our findings in the spheri- 
cal geometry, except that in our case, even at fixed en- 
ergy, the degeneracy can be lifted by the angular momen- 
tum constraint (in the microcanonical ensemble). The 
square doubly periodic domain is another example of a 
case of marginal ensemble equivalence where our analy- 
sis in terms of Goldstone mode is expected to be valid, 
with the difference that the marginal ensemble equiva- 
lence there only holds with respect to one parameter, the 
energy, whereas in the spherical case it holds in a two- 
dimensional thermodynamic parameter space, the (E, L) 
plane. 

From the geophysical application side, our study is 
based on the QG equations. These equations provide a 
simplified model of the true equations of motion for geo- 
physical fluids, which are often called "primitive equa- 
tions" and rely on the Navier-Stokes equations. In the 
QG approximation, the essential hypotheses are that 
the fluid layer is thin and that gcostrophic balance pre- 
vails approximately (which implies that rotation is strong 
enough) [37] • Both hypotheses are relatively well veri- 
fied in the atmosphere and oceans. Hence, although we 
consider here the simplest version of the QG equations 
(barotropic case) , the model in itself is relatively realistic 
(it is not a toy model). However, there is of course a fun- 
damental difference between the model and reality, as we 
consider here that the forcing and dissipation equilibrate 
locally while it is only true globally in real steady states. 
This strong simplification is common to all the statistical 
mechanics approaches up to now, as it ensures that the 
Liouville theorem is satisfied. 

It is generally believed that the dynamics of the atmo- 
sphere is essentially due to out-of-equilibrium processes. 
Nevertheless, the equilibrium states obtained here resem- 
ble to some extent to patterns observed in planetary at- 
mosphere. For instance, the solid-body rotation corre- 
sponds roughly to the first order of the general circula- 
tion. This type of flows was also used in the past to 
model the phenomenon of super-rotation in planetary 
atmospheres |31j . The similarity is much more strik- 
ing when considering the saddle points of the entropy 
functional (see [33J for a detailed description), which al- 
though not formally stable, may be long-lived as the sys- 
tem may not generate spontaneously the destabilizing 
perturbations. In particular, we find states with vortices 
in the middle latitudes which are reminiscent of the at- 
mosphere of the Earth. The lowest order saddle point is 
a quadrupole state, in agreement with simulations of 2D 
flows on a non-rotating sphere [3H1 US]- In the rotating 
case, zonal structures prevail in numerical simulations 
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[4"8l I5U] . When mixing is not strong enough, numerical 
simulations [51] and laboratory experiments [52] indicate 
a discrepancy between the final organization of the flow 
and the statistical mechanics prediction. In this case, 
the system may remain in the vicinity of saddle points of 
the entropy functional for a long time. Another example 
where the out-of-equilibrium steady-states are well de- 
scribed by the equilibrium theory [53J is that of the von 



Karman flow (the functions relating stream function, an- 
gular momentum and azimuthal velocity being selected 
by the forcing and dissipation [H]). This indicates that 
the out-of-equilibrium attractors may concentrate near 
the equilibrium states and thus computations carried out 
in the framework of the equilibrium theory may remain 
relevant to describe the out-of-equilibrium steady-states 
when forcing and dissipation equilibrate each other. 
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